Analytical solutions for a two-level system driven by a class of chirped pulses. 
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We present analytical solutions for the problem of a two-level atom driven by a class of chirped 
pulses. The solutions are given in terms of Heun functions. Using appropriate chirping parameters 
an enhancement of four-orders of magnitudes in the population transfer is obtained. 
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Interaction of coherent optical pulses with quantum 
systems is a fundamental problem [TJ [SJ. Nowadays, the 
laser systems produce controlled intense ultra-short op- 
tical pulses [3]. Various technologies have been used for 
pulse shaping [4] that allows researchers to provide co- 
herent optical control [5] of excitation in quantum sys- 
tems that has a broad range of applications from non- 
linear laser spectroscopy to generation of coherent radi- 
ation. For example, the chirped pulses [TJ |4] are used 
to produce maximal coherence in atomic and molecular 
systems. Maximal coherence can be used for generation 
of short-wavelength of radiation that has been a focus 
of research recently where an atomic system under the 
action of a far-off resonance strong pulse of laser radi- 
ation has been considered and it has been shown that 
such pulses can excite remarkable coherence on high fre- 
quency far-detuned transitions; and this coherence can 
be used for effcient generation of soft x-ray and ultra- 
violet (XUV) radiations [5J [7] . Also maximal coher- 
ence can be used for molecular spectroscopy, for exam- 
ple, time-resolved coherent Raman spectroscopy, to ob- 
tain molecule-specific signals from molecules, which can 
serve as a marker molecule for bacterial spores pQ. 

In this brief report we will investigate two class of 
chirped pulses for which the problem can be solved ex- 
actly in analytical form. Using the appropriate chirp- 
ing parameters, the population transfer, after the the 
pulse is gone, can be optimized and for the pulse consid- 
ered here, four-order of magnitudes enhancement was ob- 
tained. Unchirped pulse corresponding to Heun and Con- 
fluent Hcun equation has been recently investigated ex- 
tensively [7] where we have included an estimate of en- 
ergy of emission of soft x-ray and ultraviolet radiation via 
excited quantum coherence in the atomic system. The 
estimate shows good potential for a source of coherent 
radiation based on the discussed mechanism. 

The equation of motion for the probability amplitudes 
for the states \a) and \b) (see Fig la) of a two-level atom 
(TLA) [51 H] interacting with a classical field (under 
rotating-wave approximation RWA) with non-zero chirp- 
ing [ITj . is given as 

c a = in(ty*vc b , 



c b = in*{t) e - l§{t) c a 



(la) 
(lb) 



1 


f A 




ft 




CO 


U V n 



(a) 




|b> 



i.o I "(b) 



FIG. 1. (Color online) (a)Two-level atomic system, atomic 
transition frequency uj — ui a — cot, detuning A = u — v and 
Rabi frequency Cl(t) — p£(t)/2h. (b) Classical electromag- 
netic field E(t)= exp(-a 2 £ 2 )cos(i/i + t 2 ) 



where <&(t) = At+ <j>{t). Here A = u - v [ID] and £l(t) = 
p£{t)/2h. Let us define the dimensionless parameters as 
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To solve for C a , we can get a second order linear differ- 
ential equation for C a (t) from Eq.jfll), which in terms of 
the dimensionless parameters Eq.fpFis given as 



C a 



i/3- 



»(r) 
fi(r) 



(r) 



C a + n 2 (r)C a = 0. (3) 



Class I: Heun Equation 

To find analytical solution for Eq.©, We introduce a 
new variable tp = f(r) defined by 



= (l/2)ln[^/(l - V r +X ] 



(4) 



and make an ansatz for the pulse envelope f2(r) and the 
chirping function 4>(t) as 
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1/2 
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(c-cp) \ \fi + X(p 

-2cC + 2[(C + + c(C + ^)]v 



(5a) 
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((p-c)(n + \ip) 
In terms of the variable f(r) and the definition of 




FIG. 2. (Color online)(I) Heun Equation case: Chirping function 4>{ T ) given by Eq.(15l for £ = 10 and (a) 8 > 1, (b) 8 ~ 1 



(c) The Electric field E(r) for varying £ and 8 — 1.01. (d) Probability amplitudes for the upper level \a) for the corresponding 
fields in (c). |C"„(r)| is given by Eq.fll) 



0(r),^(r) from Eq.(|5j), Eq.(|3j) takes the form 

abip — q 
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where (c > 1) and 
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The parameters of a Heun Equation [T^l EH] are con- 
strained, by the general theory of Fuschsian equations 
as, p + a + v = a + b + 1 which provides us a the first 
constraint relation for the chirping parameters CiVit as 



C + v + C = 0. 



(8) 



The quantity ^ + </> is the instantaneous pulse frequency; 
thus </> should vanish for maximum of ^(t). From 
Eq.((5k), we get the corresponding tpo which satisfy the 
equation 

Xip 3 - (2A + p)ip 2 + c(A + 2p)ip - qu = 0. (9) 

Thus the second constraint relation for the chirping pa- 
rameters is given as 



-2cC + 2[(C + Q + c(C + >?)bo 



0. (10) 



The general solution for Eq.MJ), which has regular singu- 
larity at ip = is given in terms of Heun local solutions, 
W(ip) as, 

C„ = V^-PHlic, q + (1 - p)((c - l)<x + a + b - p + 1); 
a — p + 1, & — p 4- 1, 2 — p, a; <p] + P2 W [c, q; a, b, p, a; <p] , 

(11) 

where the constants, Vi,T > 2 can be found using the ini- 
tial conditions of the system. In the limit r — > 00, the 



population left in the level |o) can be obtained by substi- 
tuting tp — > 1 in Eq.( |ll[ ). Let us consider a simple case 
of p, = 1, A = in Eq!Q and Eq.([5k) gives 



v(t) 



1 + tanh(r) 



n(r)=7 



2y(l - <p) 

(c-tp) 



1/2 



From Eq( 12 ) the pulse takes the form 
7sech(r) 



n(r) 



i/5 — tanh(r) 



2c- 1. 



(12a) 
(12b) 

(13) 



This pulse shape serves as an excellent model for a 
smooth box pulse, by taking care of non-analyticity at 
its edges with the help of the pulse parameter 6. Using 



Eq.|9| we get <p 



± \/c 2 — c. From one of our ear- 
lier assumptions c > 1 only one of the possible values 
is allowed for ipo as < tpo < 1. Subsequently using 



ipo = c — \Jc 2 — c in Eq.(10l, we get the constraint equa- 



tion as 






= S-y/6 2 -l. 



(14) 



The defining equation for the chirping function takes the 
form 

4>(j) = {{fa - \?8 2 - l) T + ln^cosh(r) - sinh(r)]} . 

(15) 



For the pulse defined by Eq.(13), using the scaling pa- 
rameters Eq.Q and the chirping function Eq.( 15 1, Eq.([3]) 
gives 



C a (r) 



1 / 1 - 25tanhr + tanhV 

2 \ S — tanhr 



S 2 -l 



S — tanh(r) 



•i/3 



C a {r) + 



7 2 sech r 



8 — tanhr 



- it [VS 2 - 1 

C a (T) = 0. 



(16) 



Let us define the initial conditions for our system as 

C a (r ->■ -00) = 0, |C b (r->-oo)| = l. (17) 



Solution for Eq.(16), satisfying the initial conditions is 



give as 

C a = Vy l - p Hl[c, q + (1 - p)((c - 1)<t + a+ 

b — p + 1); a — p 4- 1, b — p 4- 1, 2 — p,a; ip], 



(18) 



where ip(r) is given by Eq.( 12 1) and the Heun parameters 
as . 
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The chirping parameters £, ?y and the constant V are 
given as 



C = -^(i-s + Vs 2 -i 



i(i + s-Vs 2 -i 



V = ij 



l + i(2C + 0) 



(20a) 
(20b) 
(20c) 



Here we have kept £ as a free parameter for the chirping 
function <fi(t). In Fig(2), we have considered some forms 
of the chirping function <f>(r) [see Fig(2a,2b)] for d > 1 
and (5 « 1 respectively, given by Eq.(15). Influence of 



chirping on the evolution of the probability amplitude 
for the upper level \a) in shown in Fig(2d) for the corre- 
sponding pulses in Fig(2c) 



Class II: Confluent Heun Equation 

In this section we will consider another class of pulse 
and the corresponding chirping function. Let us define 
the pulse and the chirping function as 



n(r) 



2^2^(1 -<p)^/p 

2C - 2(C + rjyp 
p + Xip 



(21a) 
(21b) 



In terms of the variable <£[t) and the definition of 
CI(t),<J>(t) from Eq.([5|, Eq.M takes the form 



d 2 C a 

dip 2 

where 

1 

u = — i 
2 



dCa (pp + q)Cg 

if — \ ) dip ip(ip — 1) 



= 0. (22) 
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The critical point which corresponds to the peak of f2(r) 
is given by 



Po = 



A + 3/i - v /(A + /i)(A + 9/x) 



2A 



(24) 



and the corresponding tq can be found using Eq.Q. At 
this point = which gives a constraint relation as 



3C(A + /i) - C\/(A + m)(A + 9m) + 



(25) 



r\ (A + 3/z - v/(A + m)(A + 9m)) = 0. 
The general solution of the Confluent Heun Equation 



Eq. ( 22 ) is given as 

C a = 7qHZ (c) [0, u-l,v-l,p,q+(l- uv)/2, p] 
+ V 2 p 1 ~ u ni < - c ' ) [0, l-u,v-l,p,q+(l- uv)/2, ip], 

(26) 

where V\, V2 can be found using the initial condition of 
the system. In the limit r — > 00, the population left in 
the level \a) can be obtained by substituting p — > 1 in 
Eq.(26). Let us consider a simple case of p = 1,A = 



0. Thus the new variable is given by Eq.(12a,). From 



Eqs.(24 25) we get p = 1/3 and r\ = 2£. Pulse shape 



Q (t) and the chirping function can be written as 

Q( T ) = 7 sech(r) [1 - tanh(r)] 1/2 , (27a) 

4,{t) = -C [1 4- 3tanh(r)] . (27b) 



For the pulse defined by Eq.pTk), using the scaling 
parameters Eq.Q and the chirping function Eq.(27:>), 
Eq.J3]| gives 

C (r)- li0 - i [1 + 3tanh(r)] - »C[1 + 3tanh(r)]| C a (r) 

+ 7 2 sech 2 (r) [1 - tanh(r)] C a {r) = 0. 

(28) 

Let us define the initial conditions for our system as 



Eq.(17). Solution for Eq.(28), satisfying the initial con 



ditions is give as 
C a = Vp 1 - u 'Hl {c) [0, 1-u, v-l,p, q+(l-uv)/2, tp], (29) 



where <p(r) is given by Eq.( 12 i) and the Heun parameters 
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FIG. 3. (Color online) (II) Confluent Heun Equation case: (a) 
Profile of the Electric field E(r) for varying f. (b) Probability 
amplitudes for the upper level \a) for the corresponding field s 
in (a). The pulse envelope f2(r) and 0(r) is given by Eq.(27l. 
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FIG. 4. (Color online) Effect of chirping on the population 
left in the excited states \a). (a) Heun Equation and (b) 
Confluent Heun Equation. The inset shows the dip (minima) 
in the population left for the Heun case. For calculations 
13 = 2.5,7 = 0.25,(5 = 1.01 (c) Chirping function for the 
Confluent Heun case d(r) 



The constant V is given as 



V = 7 



2 (3/2+3iC) 

(2C + /3)-i 



(31) 



Here we have kept £ as a free parameter for the chirping 
function </>(t) given by Eq.(27). Influence of chirping on 
the evolution of the probability amplitude for the upper 
level | a) in shown in Fig(3b) for the corresponding pulses 
in Fig(3a) . To see the effect of chirping on the population 
left in the upper level \a), we have plotted in Fig(4), 
|C a (oo)| as a function of the free chirping parameter for 
the Heun and the Confluent Heun case for a particular 
choice of the detuning /3 and the peak Rabi frequency 
7. We see that |C (oo)| ranges from 4 • 10~ 3 ^ 6 • 10 _1 . 
In Fig(5) we have plotted population left in the upper 
level I a) as a function of detuning /3 and the free chirping 
parameters £, £ for the Heun and the Confluent Heun 
case respectively. 

In conclusion, we have found new solutions for the 
problem of two-level system driven by chirped pulses in 
exact analytical form. The solutions are given in terms 
of Heun function HI and the Confluent Heun function 
Hl( c > . Using the appropriate chirping parameters we 
have shown that the population left after the pulse is 
gone can be enhanced by four-orders of magnitude [see 
Fig(4(a) inset]. Application of chirped pulses on the pop- 
ulation transfer or generation of coherence makes it very 
interesting to look for new shapes for which the system 
can be solved analytically because finding exact analyti- 
cal solutions for such a problem will not only supplement 
numerical simulations but will also be useful in under- 
standing the underlying physics. 
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FIG. 5. (Color online) Contour plot showing the effect of 
chirping and the detuning on the population left in the excited 
states I a) .(a) and (b) corresponds to Heun while (c) and (d) to 
Confluent Heun case. For calculations (a) 7 = 0.50,5 = 1.01 
(b) 7 = 1, S = 1.01 (c) 7 = 0.5 (d) 7 = 1. 
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